Multimode quantum properties of a self-imaging OPO: 
squeezed vacuum and EPR beams generation 
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We investigate the spatial quantum properties of the light emitted by a perfectly spatially de- 
generate optical parametric oscillator (self-imaging OPO). We show that this device produces local 
squeezing for areas bigger than a coherence are that depends on the crystal length and pump width. 
Furthermore, it generates local EPR beams in the far field. We show, calculating the eigenmodes of 
the system, that it is highly multimode for realistic experimental parameters. 
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I. INTRODUCTION 

Highly multiplexed quantum channels are more and 
more needed as complexity increases in the quantum 
communication and information protocols. They can be 
obtained by coupling many single mode quantum chan- 
nels [1], but also by directly using highly multimode 
quantum systems. In addition, the resolution of several 
problems in quantum imaging [2] requires the generation 
of non-classical states of light having adjustable shapes in 
the transverse plane: this is the case for superresolution 
0] , or for image processing below the standard quantum 
noise level [4]. For all these reasons, it is very important 
to develop a source of highly multimode non classical 
light (squeezed and/or entangled) of arbitrary transverse 
shape. 

In the continuous variable regime, where optical res- 
onators are necessary to efficiently produce non-classical 
states, one of the keys to successfully generate multimode 
light is the ability to operate a multimode optical res- 
onator. Indeed, many theoretical proposals rely on the 
use of an Optical Parametric Oscillator (OPO) operated 
below threshold with planar cavities [H or with confocal 
cavities [6[ [7[ [8[ which spatially filter half of the trans- 
verse modes. However, these proposals still did predict 
the arising of local vacuum squeezing and image ampli- 
fication. Hence, we propose here to keep the parametric 
process to generate non-classical light but also to over- 
come the problems encountered by the use of a full trans- 
verse degenerate cavity : the self-imaging cavity [9] . This 
type of resonator, used for instance to improve the power 
of multimode lasers [lOj], is in principle able to transmit 
any optical image within its spatial bandwidth. 

The aim of this article is to demonstrate that the self- 
imaging OPO is an excellent candidate to produce lo- 
cal squeezing, image amplification and also local EPR 



beams, taking into account its physical limitations such 
as the thickness of the crystal and the finite size of the 
various optical beams and detectors. 

The following section (section [Til) describes the experi- 
mental configuration and develops the theoretical model, 
as well as the method used to determine the squeez- 
ing spectra measured in well-defined homodyne detection 
schemes. In section HIH the results for such quantities re- 
spectively in the near field and in the far field are given, 
and we investigate the generation of EPR beams. Finally, 
in section HV1 we compute the eigenmodes of the system 
and show that they are closed to Hermite Gauss modes. 



II. SELF-IMAGING OPTICAL PARAMETRIC 
OSCILLATOR 

A. The self- imaging cavity 

We consider the parametric down conversion taking 
place in a self-imaging optical parametric oscillator whose 
cavity has been depicted in the pioneer article of Arnaud 
[9[. Such a cavity is a fully transverse degenerate one, 
which implies that all the transverse modes of same fre- 
quency resonate for the same cavity length. From a ge- 
ometrical point of view an optical cavity is self imaging 
when an arbitrary ray retraces its own path after a single 
round trip. In such a cavity, in the paraxial approxima- 
tion, the ABCD matrix M after one round trip is equal 
to identity: 



A B 
C D 



1 
1 



(1) 



The simplest self imaging ring cavity requires three lenses 
of focal length f i: i = 1, 2, 3 (@). As depicted in FigQ, 
the ring cavity is self-imaging provided the distances Cij 
of the image plane of the lens i and the objet plane of 
the lens j are given by: 
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here, the pump is not depleted, and fluctuations of the 
pump field do not contribute, at first order to the fluctu- 
ations of the signal. 

The interaction Hamiltonian of the system, taking into 
account the thickness of the crystal and the shape of the 
pump, is given by 
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where g is the coupling constant proportional to the sec- 
ond order nonlinear susceptibility - 



FIG. 1: Self-imaging OPO scheme in a ring cavity 
configuration 



C. Evolution equation of the field 



Let us consider an Optical Parametric Oscillator 
(OPO) whose cavity is the self-imaging one described in 
Fig III ([7]). A type I parametric medium of length l c is 
centered on the plane C located at the longitudinal coor- 
dinate z=0. The OPO is pumped by a gaussian TEMqq 
field E p of amplitude A p and frequency uo p = 2co s . Its 
waist w p is located at the plane C. The OPO works in 
a longitudinal degenerate operation for which signal and 
idler have the same frequency uo s . We assume that for 
the pump wave, all the mirrors are totally transparent, 
and that for the signal field, the coupling mirror M c has 
a small transmission t, the other three mirrors being per- 
fectly reflecting. 



In this section, we investigate the intracavity evolu- 
tions of the signal field in the crystal plane C (near field) 
and of its spatial fourier transform (far field), taking into 
account the crystal thickness and the finite size of the 
pump. The nonlinear interaction is supposed to be very 
weak, so that the field amplitude in a single pass through 
the crystal is only slightly affected. Therefore, the z de- 
pendence of the operators a\ can be removed in Eq©. 
The longitudinal variation of the signal operator B is 
due to the diffraction described in the modal functions 

1. Near-field evolution 



B. Electric field operators 

We will follow an operational approach close to 
the one developed in the confocal case [12], and in order 
to keep the present article concise we give only the main 
steps of the calculation. The intracavity signal field at 
frequency uo s is described by a field envelope operator 
5(x, z). In the self imaging resonator, at resonance, the 
field can be decomposed on any transverse mode basis 
(such as the Gauss-Laguerre modes for instance). The 
field operator becomes : 

B(x,z,t) = ^2fi(x.,z)ai(z,t) , (3) 
l 

where a\(z, t) is the annihilation operator of a photon in 
mode I at the cavity position z and at time t. fi is the 
amplitude of the I mode. This field obeys the standard 
equal time commutation relation at a given transverse 
plane at position z: 

[B(x,z,t),B\ X \z,t)]=5( X - X >). (4) 

Indeed, contrary to the confocal case [12[ or to any par- 
tially imaging cavity case, this operator is the same as 
the one in the vacuum as no spatial filtering is induced 
by the cavity. In the regime below threshold considered 



At the mid-point plane z = of the crystal, designated 
as the near field plane in the following, the B field evolu- 
tion can be expressed as the sum of a damping and free 
propagation term inside the cavity and of a parametric 
interaction term : 

^(x,0,t) = - 7 (l + W)B(x,0,t) (6) 

+g J J d 2 x"K int (x,x")B^(x",0,t) + V^7 s <n(x,0,t) 

where 7 is the cavity escape rate, 5 the normalized cavity 
detuning of the modes, and B{ n the input field opera- 
tor. Ki n t is the integral kernel describing the non-linear 
interaction. Assuming exact collinear phase matching 
hp — 2k s , and neglecting walk off, this kernel associates 
two points x and x" through the pump amplitude at the 
average position x + x " , and a function A(x — x") describ- 
ing the diffraction effects within the crystal. 

Xi„t(x,x") = A p (^4^)A(x-x") (7) 

with 

a/ „n ik 3 f lc ^ 2 dz' ^k s < ,,.2 , . 
A(x - x") — 4- / —r e^ 7 ' 1 (8) 
4ttI c J_ lc/2 z' 



3 



where k s = n s uo s /c is the field wavenumber, and n s the 
index of refraction at frequency uo s . It can be expressed 
in terms of the integral sine function Si(x) = J* sin ^ du 



A(x - x") = 



2ttL 



2L 



(9) 



In the thin crystal case (l c — > 0) the function A(x ± x") 
tends to the usual two-dimensional distribution <$(x±x". 

In the thick crystal case, the parametric interaction 
mixes the operators at different points of the transverse 
plane, over areas of finite extension given by the spatial 
extension of the kernel Ki nt . This extension is charac- 
terized by the width of the sine function, which define a 
coherence length : 



lcoh — 



(10) 



When |x — x"| ^> l co h, A and therefore the kernel K int 
take negligible values, there is no coupling between these 
two positions. On the other hand, when |x — x"| <C l co h 
the coupling mixes the fluctuations. Thus, we can define 
l co h as the quantum resolution of our system. 

Because of the finite size of the pump, the kernel will 
take negligible values for x + x" > w p . Therefore, we 
can define the number of transverse modes excited by the 
parametric process inside the cavity, as the ratio between 
the size of the pump, and the area defined by l co h- 



J 2 

L coh 



(ii) 



This definition relies on the classical imaging properties 
of the system. We will show in the last section of this 
article show that it is consistent with the computation of 
the eigenmodes of the system. 



2. Far-field evolution 

Let us introduce the spatial Fourier transform of the 
signal field envelope operator ([l2[) 



/ 
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B(q,z,t) = I — B(x,z,t)e-* x (12) 



Equation (J6j) becomes: 
dB 



dt 



(q,0 ) t) = -7(l + *A)B(q,0,t) 



(13) 



g J d 2 q" K int (q, q") (q", 0, t) + ^B in (q, 0, t) , 

where the coupling Kernel Ki nt (<\, q") is the Fourier 
transform of the kernel (j7J) with respect to both argu- 
ments. Straightforward calculations show that 



^mt(q,q') = A p (q + q')sinc[ 



where A p is the spatial Fourier transform of the Gaussian 

~ 2 2 

pump profile, i.e. A p (q) = ^A p exp (-|q| 2 ^). 

The sine term in the coupling kernel of Eq. ([T4|) is 
the Fourier transform of the A terms in Eq. 0), and 
correspond to the limited phase-matching bandwidth of 
the nonlinear crystal. For a thin crystal, phase match- 
ing is irrelevant and there is no limitation in the spatial 
bandwidth of down-converted modes, whereas for a thick 
crystal, the cone of parametric fluorescence has an aper- 
ture limited to a bandwidth of transverse wavevectors 
Aq « I /l co h oc 1/y/XQ. In the self imaging geometry, the 
cavity ideally transmits all the Fourier modes, so that the 
spatial bandwidth is only limited by the phase matching 
along the crystal. Finally, we have to notice that in the 
far field configuration the A p (c[ + q') couples different 
q- vectors modes within the finite width of the pump. 



D. Input /output relation 

In order to calculate the noise spectrum of the outgoing 
field, an input/output method is used. The input field is 
supposed to be in a coherent state and the fluctuations 
at the output can be inferred. The relation linking the 
outgoing fields B out (^i) to the intracavity and input 
fields at the cavity input /output portpjl is: 



B°« t (x,t) = ^i?(x,t)-B i "(x,t) 



(15) 



The evolution equation of the field, either in the near or 
in the far field can be solved in the frequency domain by 
introducing: 



dt 



B in/out^ t y 



iQt 



which lead to the input/output relation, linking 
B m (x,l]) and 5°^(x,l]). 

In the case of a thin crystal in the near field [8], this 
relation describes an infinite set of independent optical 
parametric oscillators. In this case the squeezing spec- 
trum can be calculated analytically. More generally, this 
relation in near field links all points in the transverse 
plane within the coherence area. In order to get the in- 
put/output relation, we have to inverse the input /output 
relation by using a numerical method used in jl2j ]. 



E. Homodyne detection scheme in the near field 
and far field 



In the following sections we calculate the noise spec- 
trum at the output of the OPO as a function of the de- 
tected transverse mode selected by an homodyne detec- 
tion scheme [14]. By mixing it with a coherent Local 
Oscillator (LO) of various shape on a 50% beamsplit- 
ter (reflection and transmission coefficients r = and 

t = -^), one can measure the fluctuations on any trans- 
verse mode of the output of the self imaging OPO, by 
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measuring the photocurrents difference. The two identi- 
cal detectors of different size and position are supposed to 
have a perfect quantum efficiency. All the fields are evalu- 
ated at the beam-splitter location, and the z-dependence 
is omitted in the following. 

We use two different configurations : near-field (x- 
position basis) and far-field (q- vector basis). The com- 
plete detection scheme is schematically shown in Fig. 
[2] and [3l In the near-field configuration, the imaging 
scheme is composed of a two-lens afocal system (focal 
length /) which images the crystal/cavity center plane C 
onto the detection planes D and D'(near field planes). In 
the far field configuration, a single lens of focal length / 
transforms its focal object plane C into the image focal 
detection plane D. Any image in the object plane C is 
tranformed into its fourier transform in the plane D (far 
field plane). 




LO 



FIG. 2: Balanced homodyne detection scheme in the 
near field. Two matching lenses of focal f are used to 
image the cavity center C at the detection planes D and 

D' 




FIG. 3: Balanced homodyne detection scheme in the far 
field. A matching lens of focal f is used to obtain the far 
field image of the object plane C at the detection planes 
D and D' 

For near-field imaging, the local oscillator can be ex- 
pressed as ol(x, z) = z)\e l(pL ^ z \ The difference 
photocurrent is a measure of the quadrature operator: 



E H (SI) 



=/dx[ 

J det 



B° ut ( X , OK(x) + B™ t+ (x, -fi)a L (x)|l6) 



where det is the image of the photodetection region at 
the crystal plane C, and assumed to be identical for the 
two photodetectors. The quantum efficiency of the pho- 
todetector is assumed to be equal 

For far-field imaging, the lens provides a spatial Fourier 
transform of the output field B out (x,Q), so that at the 



location of plane D the field Bjj t (x, Q) is: 
B^(x,n) = —B out (—x,Q) 



(17) 



In this plane, B^ Lt (x, Q) is mixed with an intense station- 
ary and coherent beam a^ Q (x) = j^&Loi^xjr, wnere 
oll{x) has a gaussian shape, with a waist wlo- The ho- 
modyne field has thus an expression similar to the near 
field case, where functions of x are now replaced by their 
spatial Fourier transforms: 

E H (n) = (is) 

dq[B out (ci, n)at (<i) + B° ut+ (q, -fi)a LO (q)] 

det 

In near- and far-field, the fluctuations 5Eh(Q) of the 
homodyne field around steady state are characterized by 
a noise spectrum: 

/+oo 
dSl'(SE H (Sl)SE H (Sl')) = N + S(Sl) (19) 
-co 

where Eh is normalized so that N gives the mean photon 
number measured by the detector 



N : 



det 



dx\a L (x)\" 



(20) 



N represents the shot-noise level, and S is the normally 
ordered part of the fluctuation spectrum, which accounts 
for the excess or decrease of noise with respect to the 
standard quantum level (S=0). One should note that 
there is a complete equivalence between a setup with a 
finite and flat local oscillator and infinite detectors, and 
a fiat and infinite local oscillator combine with finite size 
detectors. We will often use the configuration with finite 
size photodetectors in the following. 



III. NON-CLASSICAL PROPERTIES 

We expose here the main properties of the fields emit- 
ted by the sub-threshold self-imaging OPO. We will first 
consider the squeezing in the near-field in a very similar 
manner as what was done for a confocal OPO. Then we 
will study the far field properties and demonstrate local 
EPR correlations. 



A. Squeezing in the near field 

As the self imaging cavity does not exert any spatial fil- 
tering on the fields, the non-classical properties are very 
similar to those observed in the single pass configura- 
tion. We will here show the main squeezing predictions 
for such a device, taking into account the thickness of 
the crystal. The corresponding calculations are avalaible 
upon request to the authors. 
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Let us first consider the case of the thin crystal approx- 
imation, where no characteristic length is introduced in 
the model. We consider a thin crystal self-imaging OPO 
pumped by a gaussian beam. We look at the output 
quantum fluctuations with a pixel-like detector whose po- 
sition is varied. In figure [H the squeezing is plotted as 
a function of the detector distance from the optical axis 
for different mean powers of the pump (A p = 1 corre- 
sponding to the threshold on the axis). The squeezing is 
maximum when the detection is centered on the pump 
beam, and tends to zero far from the center. Figure [H 
shows that the squeezing increases with the total pump 
power and depends critically on its local value. Hence, 
any transverse position on the crystal acts as an indepen- 
dent OPO. 




FIG. 4: Quantum noise at zero-frequency normalized 
to the shot noise, for a pixel like detector located in the 
near field plane, as a function of the pixel distance from 
the origin s, normalized to the waist of the pump 
(s = ^-) and for different pump values 

The same behavior is observed in a configuration closer 
to actual experimental scheme using a circular detector 
whose radius Ap can be varied. Like in the previous 
case, the curves in figure [5] are crucially dependent on 
the pump power. 

In a more general realistic study, we have to take into 
account the finite size of the crystal. For a thick crystal, 
the coherence length l co h introduced in equation ([TO]) has 
to be taken into account. On transverse size smaller than 
this coherence length, fluctuations are mixed inside the 
crystal. 

In a first step, we consider a quasi plane pump for which 
its waist w pi considered as infinite, is much larger than 
the coherence area and can excite many modes. The de- 
tector is centered with respect to the optical axis and its 




FIG. 5: Quantum noise at zero- frequency, normalized 
to the shot noise, in the thin crystal case, for a circular 
detector in the near field plane, as a function of the 
radius of the detector normalized to the pump waist 
As = ^ 

UJ p 



size can be changed. As shown in figure El the squeezing 
is maximum when its size is larger than the coherence 
area, and the noise tends to shot noise for a pixel like 
detector, whose size becomes smaller than the coherence 
length. At the scale of the coherence area, the OPO can 
be considered as locally single mode and the local fluctu- 
ations are mixed. In regions smaller than the coherence 
area independent modes having their own fluctuations 
cannot be excited. For large detectors, several coherence 
area can be excited, the OPO can be considered as mul- 
timode, and the squeezing is maximum. The coherence 
length sets the limit between single mode and multimode 
operation. 

In a second step, we have to consider a more realistic 
case for which the pump waist w p is finite. In figure [71 
we represent the quantum noise as a function of the de- 
tector size for different pump waist (normalized to the 
coherence length). For detectors smaller than the coher- 
ence area, the quantum noise goes to shot noise whatever 
the size of the pump is. As explained in the last para- 
graph, at that scale the OPO can be considered as locally 
single mode and no squeezing can be obtained. For detec- 
tors whose size is close the coherence area one, squeezing 
is obtained. Nevertheless this squeezing degrades for a 
given pump waist when the size of the detector increases. 
For a given pump size, when the detector becomes larger 
than the excited surface, vacuum fluctuations are cou- 
pled to the detected signal and the squeezing degrades. 
In the same way, for a given detector size, the squeezing 
decreases with the waist of the pump. In fact a finite 
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pump size limits the number of excited modes. Increas- 
ing the pump size, increase the number of excited modes 
and improve the squeezing. 




Ap/I u 
r coh 

FIG. 6: Quantum noise at zero- frequency, normalized 
to the shot noise, in the thick crystal case, as a function 
of the detector radius (scaled to l co h).The detector in 
centered on the pump beam 



same as those in the confocal cavity, both with a plane 
pump and with a finite pump. Corresponding calcula- 
tions are also available upon request to the authors. 

The advantage of the self imaging cavity is that it does 
not couple the two symmetrical k vectors. Therefore one 
expects correlations between two symmetrical areas in 
the far field, as it will be shown in the following. 

In order to characterize the correlation level between 
symmetrical parts of the beam, we compare the quadra- 
ture field fluctuations on two symmetrical pixels. In or- 
der to get this quantities, we use the homodyne detection 
scheme (figure [8j) proposed in [14| , where two symmet- 
rical sets of two detectors measure a quadrature of the 
field at two symmetrical positions. 





FIG. 7: Quantum noise at zero- frequency, normalised to 
the shot noise, in the thick crystal case, as a function of 
the radial size of the detector scaled to l co h, plotted for 

... 2 

several values of the parameter b - 



l c oh 2 



B. Entanglement in the far field 

In the far field, the analysis has to be performed not 
in crystal plane (near field plane) but in its Fourier plane 
(far field plane). Squeezing can be observed in the far 
field when using a symmetric detector. Indeed, contrary 
to the near-field case, in the far field configuration the 
down conversion process couples two symmetric k vec- 
tors. Thus in order to recover the squeezing one needs 
a symmetric detector relative to the optical axis of the 
imaging system. The results obtained are therefore the 



FIG. 8: Homodyne detection scheme for the 
measurement of the quadrature components of the 
output field on two symmetrical pixels: pixel 1 and 
pixel 2 

Let us consider a pixel-like detector with finite detec- 
tion area Apj, according to equation fTHl the detected field 
quadrature is given by : 



4>) 



L 



B ou *(q ) n)|a L (q)|e-^- 
B°" t+ (q,-0)|a L (q)|e^ 



(21) 



where we have introduced explicitly the phase of the lo- 
cal oscillator. To compare the fluctuations of the field 
quadrature measured in two symmetrical pixels j = 1 
and j = 2, we compare the sum and the difference of 
these quantities. 



E^(Sl)=E^(Sl)±E^(Sl) 



(22) 



In order to evaluate the degree of correlation or anti- 
correlation, we introduce the corresponding fluctuations 
spectra: 

/+oo 
dn'(E£\n)E£\n>)) (23) 
-00 



Straightforward calculations show that: 



y{+) 



(SI) 



It results that the correlation between E 



(i) 



(24) 



and 



E 



(2) 



is the same that the anticorrelation between 
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the corresponding orthogonal quadrature components 
^0 1 ^ +7r /2 an ^^+7r/2- -"- n or der to calculate ([23jh we de- 
velop the expression, so as: 

/+oo 
-00 

/ + OO 
-CO 

/+00 
d0'<4>)4 2 >')> 
-co 

/ + CO 
dtf{E™(Sl)EV(tf)) (25) 
-co 

The terms f+™ d!Q! (Ef L {Q)Ef L (Vt')) correspond to the 
result of the fluctuation spectra of a homodyne detection 
scheme using a single pixel. The other terms are cross 
correlation terms, so that: 



v + 



(fi) 



(26) 



When these variances are bellow one, EPR beams are 
obtained at the output of the self-imaging cavity. One 
should note that these variances correspond to the fluctu- 
ation spectrum obtained performing a homodyne detec- 
tion in the far field with symmetric detectors : the usual 
connection between squeezing and quantum correlations 
is exhibited, both side of the same phenomenon [15[ . More 
specifically, the spatial entanglement in the far field arise 
from the correlations between the modes a q ~ e* q,x and 
a_ q ~ e _zq,x . As a q and a_ q are EPR entangled beams, 
it is well know that the combination of modes: 



V2 



V2 



cos(q.x) 
' sin(q.x) 



(27) 



will be squeezed with respect to two orthogonal quadra- 
ture components. The modes proportional to cos(q.x) 
are the even modes: using an even detection scheme it is 
possible to see squeezing, as already found in the previ- 
ous section. Note that if we use an odd detection scheme 
(symmetrical detectors with an odd local oscillator), it 
will be also possible to see squeezing in the far field, but 
on the orthogonal quadrature. 

In order to ascertain the inseparable character of this 
physical state, Duan et al.[16] have shown one needs to 
make two joint correlation measurements on non com- 
muting observables on the system. They have shown that 
in the case of gaussian states there exists a criterion of 
separability in terms of the quantity Si 2, that we will call 
'separability', and is given by: 



1 



(28) 



The suddicient Duan criterion for inseparability is given 
by: 



First, we can perform a joint correlation measurement 
using two split detectors of same but variable size as de- 
picted in Figure [9l Figure [10] shows the evolution of the 
separability at zero frequency Si2(0) = Si 2 for different 
b parameters, in function of the detector radius scaled to 
Icohf = Xf /2ttw p [12]. Notice that results are the same as 
a local squeezing measurement using a circular detector 
of variable radius Ap centered on the optical axis. 




FIG. 9: Detection scheme for the inseparability 
measurement 




Si 2 (fi) < 1 



(29) 



FIG. 10: Inseparability at zero- frequency, and at 
resonance, as a function of the radial amplitude of the 
detector Ap(scaled to the coherence area l co hf) > m the 

finite pump regime and far field approach and for 
different values of b. 



Fig. HU shows the results obtained in the case of two 
symmetrical pixels (pixel of size equal to the coherence 
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FIG. 11: Inseparability at zero- frequency, and at 
resonance, as a function of the distance between the two 
pixels p (scaled to the coherence area l co hf, i n the finite 
pump regime and far field approach and for different 
values of b. 



Number of the mode 



(b) 



FIG. 12: a) Coupling Matrix K int between two points 
of the crystal b) Spectrum of this matrix. 



IV. PIXEL-BASED MODEL FOR THE 
SELF-IMAGING OPO 

We have so far described the non-classical properties 
of the self-imaging OPO using detection based geometry, 
very appropriate to describe actual experiments. How- 
ever, it is known that any input /output system can be de- 
scribed by eigenmodes : for instance in the case of mode- 
locked pulses of light incident on a non-linear cristal, it 
has been shown that independent modes could be found, 
either using the Schmidt decomposition in the single pho- 
ton regime [17[ or diagonalising the coupling matrix in 
the continuous wave regime [l8[ • We propose here to use 
the same technique to exhibit the eigenmodes of the sys- 
tem and give a more precise value of the number of modes 
involved in the process. 

Let us pixelize the transverse space with pixels much 
smaller than l co h and develop the OPO equations onto 
the pixel operators. To simplify the system, we can first 
consider a one dimension pixelization. Let L be the size 
of the pixelized zone, and N the number of pixels. The 
pixel i is defined as the zone of size near the abscissa 
Xi = z-^r, with i ranging from —N/2 to N/2. The pixel 



operator is therefore : 

Bi = I dxB(x) (30) 

The pixel size must be chosen small enough to ensure 
the constant value of A p and Ki nt on every pixel. In this 
case, the Kernel can be written as : 

K int (i,j) = K int (xi,Xj) (31) 

and the evolution equation (|6|) at zero frequency then 
becomes : 



- ^ + ~n ? Km (i ' + ^ B * n = (32) 

3 

To solve these N coupled equations, one must find the 
eigenvectors and eigenvalues of the matrix Ki nt The 
K matrix and its spectrum are represented on figure [121 
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FIG. 13: Shape of the eigenvectors of matrix Ki nt 
(green) for the two highest positive (top) and negative 
(bottom) eigenvalues compared to the one of the pump 
(blue) 



Its diagonalization gives a set of eigenmodes with cor- 
responding eigenvalues. Some of these eigenmodes are 
represented on figure [13j they are very close to Hermite- 
Gauss polynomials shapes whose characteristic waist is 
imposed, in our case, by the pump waist. 

These modes form a basis of uncorrelated modes of 
the emitted light. Indeed, let us call Ck the eigenmode 
of eigenvalue As K int (i,j) is both self-adjoint and 
real, and Ck components are all real. In this basis, 
equation [32] can be rewritten as set of equations, one per 
mode : 



the quadratures of the eigenmode C£± £ , normalized to 
the shot-noise level, are given by : 



V k ± = A 



1 =F r 



l±r- 



(38) 



where r is the pump power normalized to the thresh- 
old and Xmax = max/c A& the highest eigenvalue of Ki nt . 
Xmax is of special interest since it is related to the pump 
power at threshold and C m ax is the corresponding las- 
ing mode. One can see in the previous equation that for 
each mode whose eigenvalue is different from zero one 
of its two variances is bellow one, implying that it is 
non-classical. However, for eigenvalues very small com- 
pared to Xmax, the squeezing is negligible. Thus one can 
compute the number of relevant mode of the system, for 
instance using a threshold eigenvalue (about 10% of the 
maximum eigenvalue). Another possibility is to calculate 
the cooperativity [17], defined from the eigenvalues of the 
matrix. 



2\2 



(39) 



The obtained number of modes is very close to the one 
defined in equation ([TT]) in a ID case. For example, using 
typical experimental values (1 cm long crystal of index 2 
and a 300 firn at 1064 nm), we find b = -p^ = 7.5 and 
k = 6.8. This means that in the 2D case, our self-imaging 
OPO can potentially excite 50 modes. 

One should note that from these eigenmodes it is pos- 
sible to find the noise properties of the pixel operators 
B° ut after the cavity, in the near field of the crystal, by 
inverting the Ki nt matrix, which gives : 



- 7 C fc + 9 -±\ k Ct + y^Cf = (33) 

These equations can again be decoupled, using the 
quadrature operators : 



Ck+ — Ck + C k 
C fe _ = -i(C k - c\) 

The final set of equations is now given by 



• 7C/C+ + ^—\ k Ck+ 



-lC k 



gL 

N 



27C& = 

2lCt = 



(34) 
(35) 



(36) 
(37) 



In this basis, using the input/output relations [15l we 
can calculate the squeezing properties of the modes Cj^, 
in the near field of the z = plane. Using the same 
method as in [18], the fluctuations at zero frequency of 



b?£ = v * c °k± = E v ^ A k c i± ( 4 °) 

(41) 



kj 



Using these expressions, we can calculate the measured 
fluctuations of the quadratures on a detector with an 
arbitrary shape : 



Vdet± 



ijEdet 



/ TDOUt TDOUt \ 

< B i± B j± > 



(42) 



These numerical simulations show the exact same results 
as the analytical results presented in section Hill 

We thus have shown two ways of solving the problem, 
each having different physical significance. Indeed, in the 
approach of section [Till we have seen that the system has 
a coherence area that defines the smallest mode having 
non-classical properties. This is relevant of quantum 
imaging applications as it gives which pixel size one 
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can address with quantum techniques. In the present 
section we have shown that a proper description of the 
system consists of an eigenmodes decomposition, modes 
that have Hermite- Gauss shape and whose squeezing 
decreases with the mode number. However, these modes 
shape are complex to measure experimentally. 
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